IN THE CLAIMS 



1. (Currently Amended) A method of inverting a 1x1 sourc e matrix, the method A n 
article comprising a machine readable medium that stores data representing a predetermined 
function, the predetermined function comprising: 

dividing the source matrix into four 2x2 sub-matrices A, B f C and D\ 

calculating a plurality of sub-matrix products from the sub-matrices; 

calculating a determinant of the source matrix dS to form a matrix determinant residue rd of 
the source matrix as rd=l/dS; 

forming a partial, inverse sub-matrix of each sub-matrix using one or more of the matrix 
products and a determinant of each sub-matrix; and 

calculating an inverse of each sub-matrix iA, iB, iC, and iD, utilizing each partial, inverse 
sub-matrix and determinant residue rd, such that an inverse of the source matrix iS is formed. 

2. (Currently Amended) The method -article of claim 1, wherein dividing the source 
matrix S into the four 2x2 sub-matrices A, B, C and D is performed according to the following rule: 



to enable storage of each sub-matrix within a pair of SIMD registers. 

3. (Currently Amended) The article method of claim 1, wherein calculating the 
plurality of sub-matrix products further comprises: 

calculating an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 

DC = adj(D)«C 

Afl = adj(A)#fl 

wherein the adj function refers to an adjoint matrix operation and the dot symbol • refers to 
a matrix multiplication operation; and 

calculating a final sub-matrix product for each of the intermediate sub-matrix products by 
computing the following equations: 




BDC = B*DC 




CAB = C •AB. 
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4. (Currently Amended) The article meteed-of claim 1, wherein calculating the matrix 
determinant residue further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD\ 
calculating a trace value by computing a following equation: 

t = trace (AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; and 
calculating a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation. 

5. (Currently Amended) The article method of claim 1, wherein forming partial-inverse 
sub-matrices further comprises: 

performing matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, BMC and 
A*dD\ and 

computing a partial inverse for each sub-matrix according to the following matrix scaling 
equations: 

pA = A* dD- BDC 
pB=C*dB- DBA 
pC=B*dC-ACD 
pD=D*dA-CAB, 

wherein pA, pB, pC, and pD reference partial, inverse sub-matrices, and the symbol * refers 
to a matrix scaling by a scalar operation. 

6. (Currently Amended) The article method of claim 1, wherein calculating an inverse 
of each sub-matrix further comprises: 

calculating an adjoint value of each partial, inverse sub-matrix pA, pB, pC, and pD, 
according to the following rules: 
iA = adj(/?A), 
= adj(/?£), 
iC = adj(/?C), 
/D = adj(/?D), 

wherein the adj() function refers to the adjoint matrix operation; 
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calculating a final sub-matrix inverse value according to the following equations: 
iA = iA * rd 
iB = iB * rd 
iC = iC * rd 
iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
forming the inverse source matrix iS according to the following rule: 



data representing a predetermined function, the predetermined function A method comprising: 
dividing a source matrix into four 2x2 sub-matrices, A, B, C and D\ 
calculating one or more intermediate sub-matrix products from one or more of the sub- 
matrices; 

calculating a determinant of the source matrix to form a determinant residue rd utilizing the 
intermediate sub-matrix products; 

scaling a determinant of each sub-matrix and the intermediate sub-matrix products using 
determinant residue rd to form final sub-matrix products; 

forming a partial inverse sub-matrix pA y pB, pC and pD for each sub-matrix using the 
scaled sub-matrix determinants and the final sub-matrix products; and 

calculating an inverse of each sub-matrix iA, iB, iC and iD, utilizing each partial inverse sub- 
matrix to form an inverse source matrix iS. 

8. (Currently Amended) The article method of claim 7, wherein calculating the matrix 
determinant residue further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD\ 
calculating a trace value by computing a following equation: 

t = trace (AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; 
calculating a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
calculating the determinant residue rd according to the following rule: 

rd = 1/dS. 



iS = 




7. 



(Currently Amended) An article comprising a machine readable medium that stores 
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9. (Currently Amended) The article methed-of claim 7, wherein scaling by the 
determinant residue further comprises: 

multiplying each determinant by the determinant residue rd according to the following rules: 
dA = dA* rd 
dB = dB * rd 
dC = dC * rd 
dD = dD * rd; 

multiplying each intermediate sub-matrix product AB and DC by the determinant residue 
rd, according to the following equations: 

DC = DC*rd 

AB = AB* rd; and 

calculating a final sub-matrix product for each of the intermediate matrix products by 
computing the following equations: 

BDC = B*DC 

DBA = D»adj(AB) 

ACD = A«adj(l)c) 

CAB = C •AB. 

10. (Currently Amended) The article method of claim 7, wherein calculating an inverse 
of each sub-matrix further comprises: 

generating an adjoint of each partial, inverse sub-matrix by computing the following 
equations: 

iA = adj(pA) 

iB = adj(pB) 

*C = adj(pC) 

iD = adj(p£>); and 

forming the inverse source matrix iS according to the following rule: 

fiA iB\ 
iS = ^ I 
{iC iDJ . 

1 1 . (Original) A computer readable storage medium including program instructions that 
direct a computer to function in a specified manner when executed by a processor, the program 
instructions comprising: 

dividing the source matrix into four 2x2 sub-matrices A, B, C and D; 
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calculating a plurality of sub-matrix products from the sub-matrices; 
calculating a determinant of the source matrix dS to form a matrix determinant residue rd of 
the source matrix as rd=l/dS; 



products and a determinant of each sub-matrix; and 

calculating an inverse of each sub-matrix iA, iB, iC, and /D, utilizing each partial, inverse 
sub-matrix and determinant residue rd, such that an inverse of the source matrix iS is formed. 

12. (Original) The computer readable storage medium of claim 11, wherein dividing the 
source matrix S into the four 2x2 sub-matrices A, B, C and D is performed according to the 
following rule: 



to enable storage of each sub-matrix within a pair of SIMD registers. 

13. (Original) The computer readable storage medium of claim 11, wherein calculating 
the plurality of sub-matrix products further comprises: 

calculating an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 



A5 = adj(A)»5 

wherein the adj() function refers to an adjoint matrix operation and the dot symbol • refers 
to a matrix multiplication operation; and 

calculating a final sub-matrix product for each of the intermediate sub-matrix products by 
computing the following equations: 

BDC = B*DC 



CAB = C •AB. 

14. (Original) The computer readable storage medium of claim 11, wherein calculating 
the matrix determinant residue further comprises: 

computing a determinant of each sub-matrix dA } dB, dC and dD\ 
calculating a trace value by computing a following equation: 



forming a partial, inverse sub-matrix of each sub-matrix using one or more of the matrix 
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t = trace (AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; and 
calculating a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation. 

15. (Original) The computer readable storage medium of claim 11, wherein forming 
partial-inverse sub-matrices further comprises: 

performing matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, B*dC and 
A*dD; and 

computing a partial inverse for each sub-matrix according to the following matrix scaling 
equations: 

P A = A*dD-BDC 
P B = C*dB-DBA 
P C = B*dC~ACD 
P D = D*dA-CAB, 

wherein pA y pB, pC, and pD reference partial, inverse sub-matrices, and the symbol * refers 
to a matrix scaling by a scalar operation. 

16. (Original) The computer readable storage medium of claim 11, wherein calculating 
an inverse of each sub-matrix further comprises: 

calculating an adjoint value of each partial, inverse sub-matrix pA y pB, pC, and pD, 
according to the following rules: 
iA = adj(pA), 
/5 = adj(/?5), 
ZC = adj(pC), 
/D = adj(pD), 

wherein the adj() function refers to the adjoint matrix operation; 

calculating a final sub-matrix inverse value according to the following equations: 

iA = iA * rd 

iB = iB * rd 

iC = iC * rd 

iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
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forming the inverse source matrix iS according to the following rule: 
OA iB\ 



iS = 



iC iD 



17. (Original) The computer readable storage medium including program instructions 
that direct a computer to function in a specified manner when executed by a processor, the program 
instructions comprising: 

dividing a source matrix into four 2x2 sub-matrices, A, B f C and D; 
calculating one or more intermediate sub-matrix products from one or more of the sub- 
matrices; 

calculating a determinant of the source matrix dS to form a determinant residue rd of the 
source matrix utilizing the intermediate sub-matrix products and the sub-matrix determinants; 

scaling a determinant of each sub-matrix and the intermediate sub-matrix products using 
determinant residue rd to form final sub-matrix products; 

forming a partial inverse sub-matrix pA, pB, pC and pD for each sub-matrix using the 
scaled sub-matrix determinants and the final sub-matrix products; and 

calculating an inverse of each sub-matrix /A, iB, iC and iD, utilizing each partial inverse sub- 
matrix to form an inverse source matrix iS. 

18. (Original) The computer readable storage medium of claim 17, wherein calculating 
the matrix determinant residue further comprises: 

computing a determinant of each sub-matrix dA, dB, dC and dD\ 
calculating a trace value by computing a following equation: 

t = trace (AB • DC) ; 
wherein a dot symbol • refers to a matrix multiplication operation; 
calculating a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
calculating the determinant residue rd according to the following rule: 

rd = 1/dS. 
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19. (Original) The computer readable storage medium of claim 17, wherein scaling by 
the determinant residue further comprises: 

multiplying each determinant by the determinant residue rd according to the following rules: 
dA = dA * rd 
dB = dB* rd 
dC = dC * rd 
dD = dD * rd ; 

multiplying each intermediate sub-matrix product by the determinant residue rd, according 
to the following equations: 

DC=DC*rd 

AB = AB*rd; and 

calculating a final sub-matrix product for each of the intermediate matrix products by 
computing the following equations: 

BDC = B*DC 

DBA = D»adj(A£) 

ACD = A»adj(/3c) 

CAB = C »AB. 

20. (Original) The computer readable storage medium of claim 17, wherein calculating 
an inverse of each sub-matrix further comprises: 

generating an adjoint of each partial, inverse sub-matrix by computing the following 
equations: 

iA = adj(/?A) 
iB = adj(pfl) 
iC = adj(pC) 
iD = adj(pZ)); and 
forming the inverse source matrix iS according to the following rule: 
(iA iB\ 
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21. (Currently Amended) An apparatus, comprising: 
a processor having circuitry to execute instructions; 

a plurality of SIMD data storage devices coupled to the processor, the SIMD data storage 
registers to store p airs of floating point vectors during matrix calculation; 

a storage device coupled to the processor, having sequences of instructions stored therein, 
which when executed by the processor cause the processor to: 

divide the source matrix into four 2x2 sub-matrices A, B, C and D\ 

calculate a plurality of sub-matrix products from the sub-matrices; 

calculate a determinant of the source matrix dS to form a determinant residue rd of 
the source matrix as rd=l/dS\ 

form a partial, inverse sub-matrix of each sub-matrix using one or more of the 
matrix products and the determinant of each sub-matrix; and 

calculate an inverse of each sub-matrix iA, iB, iC, and iD, utilizing each partial, 
inverse sub-matrix and determinant residue rd, such that an inverse of the source matrix iS is 
formed. 

22. (Original) The apparatus of claim 21, wherein the instruction to calculate the 
plurality of sub-matrix products further causes the processor to: 

calculate an intermediate sub-matrix product for each sub-matrix by computing the 
following matrix equations: 

Z>C = adj(Z>)#C 

AB = adj(A)#fl 

wherein the adj() function refers to an adjoint matrix operation and the dot symbol • refers 
to a matrix multiplication operation; and 

calculate a final sub-matrix product for each of the intermediate sub-matrix products by 
computing the following equations: 

BDC = B^DC 

DBA = D#adj(A#) 

ACD = A*adj(z>c) 
CAB = C •AB. 

23. (Original) The apparatus of claim 21, wherein the instruction to calculate the matrix 
determinant residue further causes the processor to: 

compute a determinant of each sub-matrix dA y dB, dC and dD\ 
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calculate a trace value by computing a following equation: 

t = trace (AS* DC); 
wherein a dot symbol • refers to a matrix multiplication operation; and 
calculate a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation. 

24. (Original) The apparatus of claim 21, wherein the instruction to perform matrix 
scaling further causes the processor to: 

perform matrix scaling of a determinant of each sub-matrix as D*dA, C*dB, B*dC and 
A*DdD\ 

compute a partial inverse for each sub-matrix according to the following matrix scaling 
equations: 

pA = A*dD-BDC 
pB = C*dB- DBA 
pC=B*dC-ACD 
P D = D*dA-CAB, 

wherein pA, pB, pC, and pD reference partial, inverse sub-matrices and the symbol * refers 
to a matrix scaling by a scalar operation. 

25. (Original) The apparatus of claim 21, wherein the instruction to calculate an inverse 
of each sub-matrix further causes the processor to: 

calculate an adjoint value of each partial, inverse sub-matrix pA, pB, pC, and pD, according 
to the following rules: 

iA = adj(/?A), 
iB = adj(pB), 
/C = adj(pC), 
/D = adj(/?£>), 

wherein the adj() function refers to the adjoint matrix operation; 
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calculate a final sub-matrix inverse value according to the following equations: 
iA = iA * r d 
iB = iB * rd 
iC = iC * rd 
iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
form the inverse source matrix iS according to the following rule: 

(iA iB\ 
iS = _ _ 
\iC iDj . 

26. (Currently Amended) An apparatus, comprising: 
a processor having circuitry to execute instructions; 

a plurality of SIMD data storage devices coupled to the processor, the SIMD data storage 
registers to pairs of floating point vectors during matrix calculation; 

a storage device coupled to the processor, having sequences of instructions stored therein, 
which when executed by the processor cause the processor to: 

divide a source matrix into four 2x2 sub-matrices, A, B, C and D\ 
calculate one or more intermediate sub-matrix products from each of the sub- 
matrices, 

calculate a source matrix dS to v form a determinant residue rd utilizing the 
intermediate sub-matrix products, 

scale a determinant of each sub-matrix and the intermediate sub-matrix products 
using determinant residue rd to form final sub-matrix products, 

form a partial inverse sub-matrix pA, pB, pC and pD for each sub-matrix using the 
scaled sub-matrix determinants and the final sub-matrix products, and 

calculate an inverse of each sub-matrix iA, iB, iC and iD, utilizing each partial 
inverse sub-matrix to form an inverse source matrix iS. 

27. (Currently Amended) The apparatus s ystemof claim 26, wherein the instruction to 
calculate the source matrix determinant residue further causes the processor to: 

compute a determinant of each sub-matrix dA, dB, dC and dD\ 
calculate a trace value by computing a following equation: 

t = trace (AB • DC) 
wherein a dot symbol • refers to a matrix multiplication operation; 
calculate a determinant of the source matrix dS by computing a following equation: 

dS = dA*dD + dB*dC - t 
wherein the symbol * refers to a scalar multiplication operation; and 
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calculate the determinant residue rd according to the following rule: 
rd = 1/dS. 



28. (Currently Amended) The system apparatus of claim 26, wherein the instruction to 
scale by the determinant residue further causes the processor to: 

multiply each determinant by the determinant residue rd according to the following rules: 

dA = dA * rd 

dB = dB * rd 

dC = dC * rd 

dD = dD * rd) 

multiply each intermediate sub-matrix product AB and DC by the determinant residue rd, 
according to the following equations: 

DC = DC*rd 

AB = AB*rd; and 

calculate a final sub-matrix product for each of the intermediate matrix products by 
computing the following equations: 



CAB = C •AB. 

29. (Currently Amended) The system apparatus of claim 26, wherein the instruction to 
calculate an inverse of each sub-matrix further causes the processor to: 

generate an adjoint of each partial, inverse sub-matrix by computing the following 



iA = adj(pA) 

/5 = adj(p5) 

ZC = adj(pC) 

/D = adj(/?D); and 

form the inverse source matrix iS according to the following rule: 

(iA iB\ 
iS = _ _ 
\iC iD) . 
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BDC 



B^DC 




equations: 



Please add the following new claims: 

-30. (New) A method comprising: 

dividing a source matrix into four 2x2 sub-matrices A, B, C and D; 
storing each two element row of each 2x2 sub-matrix within a single instruction multiple 
data (SIMD) register; 

forming a partial, inverse sub-matrix of each sub-matrix using one or more of a plurality of 
sub-matrix products calculated from the sub-matrices and a determinant of each sub-matrix within 
one or more SIMD registers; and 

calculating an inverse of each sub-matrix iA, iB, iC and iD, utilizing each partial, inverse sub- 
matrix and a determinant residue rd calculated from the source matrix, such that an inverse of the 
source matrix iS is formed within the one or more SIMD registers. 

31. (New) The method of claim 30, wherein forming the partial inverse sub-matrix 
further comprises: 

calculating the plurality of sub-matrix products from the sub-matrices; and 
calculating the determinant of the source matrix Ds to form the matrix determinant residue 
rd of the source matrix as rd=l/Ds. 

32. (New) The method of claim 30, wherein dividing the source matrix S into the four 
2x2 sub-matrices A, 5, C and D is performed according to the following rule: 

(A B\ 



to enable storage of each sub-matrix within a pair of SIMD registers. 

33. (New) The method of claim 31, wherein calculating an inverse of each sub-matrix 
further comprises: 

calculating an adjoint value of each partial, inverse sub-matrix pA, pB, pC, and pD, 
according to the following rules: 
iA = adj(pA), 
iB = adj(pfi), 
iC = adj(pC), 
iD = zd](pD), 

wherein the adj() function refers to the adjoint matrix operation; 
42390P13115 16 10/038,395 




calculating a final sub-matrix inverse value according to the following equations: 
iA = iA * rd 
iB = iB * rd 
iC = iC * rd 
iD = iD * rd, 

wherein the symbol * refers to a matrix scaling by a scalar operation; and 
forming the inverse source matrix iS according to the following rule: 
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